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■ Abstract 
O ■ 

We consider the semiclassical limit of quantum systems with a Hamiltonian given by the 
Weyl quantization of an operator valued symbol. Systems composed of slow and fast degrees 
of freedom are of this form. Typically a small dimensionless parameter e <C 1 controls the 
separation of time scales and the limit e — > corresponds to an adiabatic limit, in which the slow 
and fast degrees of freedom decouple. At the same time e — > is the semiclassical limit for the 
slow degrees of freedom. In this paper we show that the e-dependent classical flow for the slow 
degrees of freedom first discovered by Littlejohn and Flynn LF91 , coming from an e-dependent 
classical Hamilton function and an e-dependent symplectic form, has a concrete mathematical 
and physical meaning: Based on this flow we prove a formula for equilibrium expectations, an 
Egorov theorem and transport of Wigner functions, thereby approximating properties of the 
quantum system up to errors of order e 2 . In the context of Bloch electrons formal use of this 
classical system has triggered considerable progress in solid state physics X CNl O]. Hence we 
discuss in some detail the application of the general results to the Hofstadter model, which 
describes a two-dimensional gas of non-interacting electrons in a constant magnetic field in the 
tight-binding approximation. 
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1 Introduction 

Semiclassical approximations play an important role in the understanding of many quantum me- 
chanical phenomena. The standard situation is a Hamiltonian H = H(x, — ihVx), acting on L 2 (M"), 
that is the Weyl quantization of a real valued function H(q,p) on the classical phase space R 2 ™. 
The goal of semiclassical methods is to approximate properties related to the quantum mechani- 
cal operator H on L 2 (R") in the limit h — > using the classical Hamiltonian H on R 2 ™ and its 
flow. Examples of such properties are the spectrum and eigenfunctions of H , special solutions of 
the time-dependent Schrodinger equation or the full unitary group e~ lH i or statistical expectation 
values tr(f(H)a). 

In many situations only some physical degrees of freedom behave semiclassically, the simplest 
situation being particles with spin. Then the quantum mechanical state space is L 2 (R™, C L ), or more 
general, L 2 (M. n ,H{) :— H, where T-Li is the state space of the "fast" or "fibre" degrees of freedom. 
Then the Hamiltonian is often of the form H = H(x,—ieV x ), where now H : R 2 ™ — > £(Hf) is 
a function on classical phase space taking values in the linear self-adjoint operators on Jif. The 
physical meaning of the small parameter e < 1 depends on the concrete problem. Also in this 
setting semiclassical approximations turned out very successful and a vast literature exists. The 
main content of this work are two results which, to our knowledge, have been neither proved nor 
claimed in this generality, but have been used already quite successfully in concrete problems, in 
particular in solid state physics, see jXCNIOj and references therein. 

We now explain our main results, postponing technical details to Section 2 and beyond. The first 
question is, which classical Hamiltonian system should be related to an operator-valued function 
H(q,p) on classical phase space? We give the following answer: to each isolated non-degenerate 
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eigenvalue eo(q,p) of H(q,p) depending smoothly on (q,p) we associate an almost invariant subspace 
n e L 2 (R™, TL{ ), IF an orthogonal projection depending on e, and a specific classical Hamiltonian 
system. The restriction of the Hamiltonian H to the range of IF is a generalized adiabatic approx- 
imation which has been studied in great detail, see e.g. [Teu03] for an overview. The restricted 
Hamiltonian n e _HTF is then amenable to semiclassical approximation using the mentioned classical 
Hamiltonian system. At leading order the corresponding scalar Hamiltonian h : R 2 ™ — > R is just 
given by the eigenvalue itself, h(q,p) — eo(q,p), and this is basically well known. However, in many 
applications next to leading order effects play an important role and we show how to incorporate 
the first order corrections into a modified e-dependent classical system. The corrected Hamiltonian 
function takes the form 

h = e + eM, 

where the function M is given in terms of the spectral projection ira(q,p) of H(q,p) corresponding 
to the eigenvalue eo(q,p). In addition one needs to modify also the symplectic form on R 2 ", 

tu e := loq + e fl , 

where ujo is the standard symplectic form on R 2 ". With z — (q,p) and a, (3 — l,...,2n, the 
correction term fi takes the form 

Map := -itrH f (7To[<9 ZQ 7ro,<9 2(3 7r ]) . 

As we will explain, is the curvature form of the induced connection, often called the Berry 
connection, on the line bundle over R 2 ™ defined by the rank one projections 7To(q,p). 

The Hamiltonian equations of motion uj E (Xh,-) — dh imply that the Hamiltonian vector-field 
Xh is given by X% = —{uj e ) a Pdph and the classical equations of motion thus have the form 



q \ _ ( -eWP E n + eOP« \ / d q h 
p J ~ \ -E n + eQ qp -eQ.il ) I d p h 



+ 0(e 2 ). (1) 



Here E n = diag(l, . . . , 1) denotes the n x n unit matrix and the error term results from the fact 
that the matrix (w £ ) Q ^ was inverted only approximately. We denote the corresponding classical 
flow by (f>l and the Liouville measure induced by the symplectic form oj e by X £ . 

We prove the following two statements about semiclassical approximations in terms of the clas- 
sical system defined above. 

Stationary expectations (cf. Theorem [T]) 

Let / : R — >• IR be real valued and a : R 2 ™ — > R be a suitable scalar observable. Then 

tv H = ^_L_ (7 dA e f(h(q,p))a(q,p) + 0(e 2 )) . (2) 

Here one could think of / being a distribution function like the Fermi-Dirac distribution and a some 
trace-class observable. A nice application of this formula where the leading order approximation 

tr H (ff f(H) a) = -1— ( [ dqdp f(e (q,p)) a(q,p) + 0(e V 

would give the wrong answer, namely zero, is the derivation of a formula for the orbital magneti- 
zation in periodic media, cf. |GA03[ IXCNIOj and Corollary [3] in Section[7] 
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Egorov theorem (cf. Theorem [5]) 

One can approximate the time evolution of quantum mechanical observables using the classical flow 
<fil. Let again a : M 2n ->1 be a suitable scalar observable, 



its quantum mechanical Heisenberg time-evolution and a(t) :— a o 0* its classical evolution. Then 



Let us try to give a heuristic meaning to the modifications in the classical Hamiltonian system. 
As we will see, the projection IF is an order e modification of tt . Thus the energy of the fast 
degrees of freedom for states in the range of IP is not exactly eo, but must be modified by eM. 
This interpretation is very natural given the symmetric form of M in (|19[) . The correction to the 
geometry of phase space, expressed in the modified symplectic form uj e is necessary, since the wave 
functions in the range of IF — ttq + 0(e) take values, micro-locally, in a line bundle with nontrivial 
geometry defined by ttq. 

The literature on semiclassical approximations for Hamiltonians with operator-valued symbols 
is huge and we give a short survey in Section [5] At this point let us mention what we believe are 
the most important contributions. 

The algebraic part of the construction of the projection IF is a central step that was done by 
different approaches and independently by Helffcr and Sjostrand HS90 j and Emmrich and Weinstein 
EW96J. Emmrich and Weinstein also showed that h appears in the transport equation for the 
amplitude of the WKB approximation, but they did not make use of the Hamiltonian flow of h. 
However, the M term and the modified symplectic form ui e appeared already earlier in a work of 
Littlejohn and Flynn [LF91 , who used the equations of motion ([1]) in order to construct Lagrangian 
submanifolds on which WKB wave functions are supported. But they arrive at these equations along 
a slightly different route and, as far as we understand their work, they considered them rather a 
technical tool than the correct physical equations for the slow degrees of freedom. 

As the equations of motion for the position and quasi-momentum of a Bloch electron, a special 
case of {1} appeared for the first time explicitly in a work of Chang and Niu |CN96] . Since then 
the "modified semiclassical model" was at the basis of a lot of theoretical progress in solid state 
physics, see [XCNIOj for a review. The only rigorous justification appeared in the case of Bloch 
electrons, where an Egorov theorem in the above explained sense was proven in |PST03j and slightly 
generalized in [DLllj . 

After the success of the modified semiclassical model in solid state physics the literature con- 
cerned with its generalizations and justifications became huge. However, to our best knowledge the 
claims ([2|) and (|3|) have not been clearly stated, nor rigorously or even just systematically derived 
until now, the one exception being [PST03] . In addition to being new, we believe that our derivation 
of ^ and © is particularly simple and transparent, granted that one is willing to work with the 
Weyl calculus. 

A key feature that distinguishes our approach and |EW96] from the majority of other approaches, 
starting with |LF91j and including jPST03j . is that we make no use of eigenfunctions ipo(q,p) of 



A(t) := e : 





(3) 



and, as a consequence, for any state p = HpH that is trace-class, 



tv n (pA{t)) = tr n (pa(t)) + 0{e 2 ). 
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H(q,p), but only of the projections n (q,p). While the former are gauge-dependent, since with 
PailiP) a l so Q 1 ^ q,v ' ( Po{ ( liP) f° r real- valued / is an eigenfunction, the projections are unique. More- 
over, in some applications like magnetic Bloch electrons, a global smooth choice for <po(q,p) simply 
does not exist for geometric reasons. Exactly this problem was our motivation for developing the 
new approach presented here and we discuss the application to magnetic Bloch bands of the Hof- 
stadter model in Section [7J The general situation of magnetic Bloch electrons is the content of a 
future paper |ST12j . 

Let us stress that in this work we make no serious attempt to achieve greatest possible generality. 
Instead we focus mostly on structural aspects and on covering the application to the Hofstadter 
model discussed in Section [7] Therefore we will make stronger assumptions than presumably nec- 
essary in order to avoid distracting technicalities. But we will comment in the following on where 
and how one can relax the assumptions. Most works, including |LF91[ |EW96| . are concerned with 
asymptotic expansions of formal symbols only and ignore the problem of turning the algebraic 
computation into statements about norm convergent approximations. Note that these algebraic 
computations and thus the structure of the classical Hamiltonian system remain exactly the same 
in all applications. Our restrictive assumptions come only into play in the steps where symbolic 
expansions are turned into statements about operators and traces. If one is willing to forego 
mathematical rigor in this step, one can freely apply our results even in cases where the precise 
mathematical assumptions are not satisfied. For an alternative but much more involved proof of 
the Egorov theorem with a slightly larger range of validity we refer to [Still) . 

We end the introduction with a short overview of the structure of this paper. In Section[3]we first 
explain the precise setting and assumptions and give a brief introduction to Weyl calculus. Then 
we state the known result concerning the existence and properties of the projection II 6 defining 
the almost invariant subspace related to the eigenvalue eo- In Section [3] we introduce the classical 
Hamiltonian h and show how it is related to the restriction of H to the range of H £ . Section [4] 
contains the mathematical statements and proofs of our main results ([2]) and <j3j) . In Section [5] 
we prove that u e is indeed a symplectic form and compute an explicit formula for the Liouville 
measure X E . Then we show how Q is related to the curvature of the Berry connection. In Section [6] 
we discuss some related literature in view of our results. Finally, in Section [71 we discuss a tight- 
binding model for a two-dimensional gas of non-interacting particles in a constant magnetic field, 
the so called Hofstadter model. Applying our abstract results, we compute the free energy per 
volume at fixed temperature and chemical potential and its derivative with respect to the magnetic 
field, the magnetization. Moreover, we also compute the Hall current at zero temperature when the 
Fermi level is in a gap, leading to the quantized conductivity in terms of Chern numbers. Although 
similar results have been established by other techniques, see |SBT12j . we believe that the derivation 
of these expressions in terms of the semiclassical model is very simple and transparent, once the 
model is established. 

Acknowledgements. We thank Giuseppe De Nittis, Omri Gat, Max Lein, Gianluca Panati and 
Hermann Schulz-Baldes for inspiring discussions and Yuri Kordyukov for pointing out several typos 
and glitches in an earlier version. S.T. is also grateful to Giuseppe De Nittis and Gianluca Panati 
for sending us their unpublished notes on the Harper and the Hofstadter model. Finally S.T. thanks 
Peter Vastag for his involvement in a joint project on the semiclassics of the Hofstadter model. 
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2 The adiabatic approximation 



2.1 General setting and precise assumptions 

We assume that the state space H decomposes as H = L 2 (R n ) ® Hi, where the subscript f stands 
for "fast" degrees of freedom or for "fibre" . Using Lebesgue-Bochner spaces we may rewrite this as 

H = L 2 {R n )®H f = L 2 {R n ,Hi). 

One can replace K™ by a flat n-dimensional torus with only very minor changes to the following 
analysis. Most importantly, one has to use an appropriately modified pseudo-differential calculus 
as explained, e.g., in |PST03| . We will give some more details on the modifications when we discuss 
the Hofstadtcr model in Section [7] For manifolds with non-vanishing curvature the first order 
corrections cannot be implemented using a modified flow, see |Lam09| . 

Assumption 1. On T-L we study a self-adjoint e-dependent Hamiltonian H that is given as the 
e-Weyl quantization H = op 1 ^ (H) = "H(e, x, — ieVz) " of a semiclassical symbol 

H{z) = H (z) + eHi(z) + £ ■ z with H , H 1 e S°(B(H { )) and ^eM 2 ". 

For any Banach space X and k € K let S k (X) denote the smooth functions from M 2 ™ to X that 
are bounded together with all their derivatives by a constant times (z) := (1 + |-z| 2 ) 7 , i.e. for 
f£S k (X) 

||/|| fer := max sup || (z) k d%f(z)\\x < oo for all r G N . 
\a\ < r 

The family of seminorms || • ||^ r turns S (X) into a Frechet space. The space of uniformly bounded 
functions A : [0,e ) -> S k (X) is denoted by S k (e,X). 

Note that compared to the introduction we allow for a subprincipal symbol Hi and explicitly 
separate a linear part in the symbol for technical reasons. As will be recalled in the next section, 
Hq and Hi are bounded operators. As a consequence H is self-adjoint on the maximal domain 
of the operator £ 9 • x — ie£ p • V x denoted by D(H) in the following. The assumption that the 
symbol H is bounded up to a part linear in z, allows to prove quite directly uniform statements 
in the following. However, all results are local in phase space R 2ra and more general Hamiltonians, 
e.g. with quadratically growing symbols, can be brought, at least formally, into the required form 
by appropriate micro-localization in phase space. A rigorous implementation of this idea in the 
general case is presumably quite technical, but has been done successfully e.g. in the context of the 
Born-Oppenhcimcr approximation, see |MS09i |TW12| . For the construction of n e it is also known 
that one can allow for H to take its values in the unbounded operators on Hf, c.f. PST03 . We 
show in [ST 12] that also the other steps of our semiclassical approximation work in this case. 

Assumption 2. We assume that Ho(q,p) has a non-degenerate eigenvalue eo(q,p) such that eo : 
R 2 " — > R is continuous and satisfies the uniform gap condition 

dist(e (g,p), a(H Q (q,p)) \ {e (q,p)}) > g > . 

We denote the eigenprojection of H Q (q,p) to the eigenvalue e (q,p) by iro(q,p), i.e. 

H {q > p)TT (q,p) = e Q (q,p)TT (q,p) ■ 
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The gap condition implies smoothness of eo and of ttq and is also crucial for the construction 
of super-adiabatic subspaces, cf. Proposition [TJ It is known that at eigenvalue crossings even a 
leading order approximation in terms of a classical flow becomes much more intricate, see e.g. 
[FKG021 ILT051 IDFJ09] and references therein. 

2.2 Key formulas from Weyl calculus 

In this section we collect for later reference some key formulas and results of the Weyl calculus and 
introduce our notation. Readers familiar with the Weyl calculus can skip this section. For more 
details on e-pseudo-differential operators we refer e.g. to |Mar02, DS99 . For a short summary on 
Weyl calculus with operator- valued symbols the readers can also consult Appendix A of |Teu03] . 

For a Schwartz function ip £ S(W l ,7i[) and A 6 S (e, B(Hf)) the action of the Weyl quantization 
A = op w (A) of A can be defined through the oscillatory integral 

(AiP)(x) = — L- / dpdyA{e,\{x + y), P )^-yy^{y). (4) 

For k > 0, by the Calderon-Vaillancourt theorem, A can be extended to a bounded operator on 
H = L 2 (R ra , Hf ). More precisely there exists a constant c n depending only on the dimension n, 
such that for all e £ [0, £q) 

\\Mb(-h) < C n \\A(e)\\ ,2n+l ■ 

The composition of operators induces a composition of symbols. For any A £ S kl (e, B(Hf)) and 
B £ S k * (e,B(Hf)), there exists a symbol C £ S kl+k2 (e, B(H { )) denoted by C = A#B such that 
AB — C. The bilinear map # : S kl x S k ' 2 — > S kl+k2 is called the Moyal product and it is continuous 
with respect to the Frechet topologies uniformly in e, i.e. for any r £ No there is a f £ No and a 
constant c r ,f < oo such that 

||(A#B)( e )|| fcl+teir < c r \\A(s)\\ kl M\B(e)\\ k ^ 

for all e £ [0, £o)- The last statement follows e.g. from inspecting the proof of Thm. 2.41 in [Fol89 . 
If for A £ S k (e,X) there exists a sequence (Aj)j^ in S k (e, X) such that 



sup e- {m+1) [A{e) -J^^A^e] 



< OO 

fe.r 



for all r e N and m £ N , then one writes A x Y^jLo £3 A j in S k (e,X). If A £ S k (e, X) has an 
asymptotic expansion with coefficients Aj £ S k (X) not depending on e, then A is called a classical 
symbol, A$ its principal symbol and A\ its subprincipal symbol. 

The Moyal product C := A#B £ S kl+k2 (e, B(H t )) of symbols A £ S kl {s,B{U{)) and B £ 
S k2 (e,B(T-L()) has an explicit asymptotic expansion C X Y^jLo^Q m S kl+k2 (e, X) such that the 
remainder maps 

m 

R m+1 : x S k2 -> S kl+k2 , (A, B) ^ R m+l := £-( m+1 ) (c(e) - ^ e?Cj{e)) 

3=0 
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are continuous. The expansion starts with the pointwise product Cq(e) = A(e)B(e) and the Poisson 
bracket C\{e) — — ^{A(e), B(e)}, where 

n 

{A,B} := d p A ■ d q B - d q A ■ d p B := £ (d P] Ad qj B - d q] Ad P] B) . 

i=i 

For classical symbols A and B the Moyal product C := Aj^B is also a classical symbol with an 
asymptotic expansion starting with 

A#B x A B + e(A 1 B + A B 1 -^{A ,B }) + 0{e 2 ) . 

Since A and B are operator valued functions, they do not commute in general and neither do their 
derivatives. Hence, in general, {A, A} ^ 0, but, if the derivatives of A are trace-class, 

tr Wf ({A,A}) = 0, (5) 

because of the cyclicity of the trace. For later reference we state also the resulting formula for triple 
products of classical symbols, 

x A B C Q + eA^BqCq + eAqB^o + eAoBod (6) 
- f {A {B , C } + {A , B Q }Co + {A \B Q \C }) + 0(e 2 ) . 

Here and in the following we use the shorthand 

{A \B \C } := d p A ■ B d q C - d q A ■ B d p C . (7) 

If A — alfi { is a scalar multiple of the identity, then A and all its derivatives commute with any B. 
As a consequence one can show that in this case 

A #B - B #A x -ie{A Q , B } + 0(e 3 ) . (8) 

The fact that the remainder term in ((8]) is of order e 3 and not only e 2 is at the basis of our higher 
order semiclassical approximations. It distinguishes the Weyl quantization from other quantization 
rules. 

We will also be interested in taking traces of pseudo-differential operators and of their symbols. 
Denoting by J\(Hi) C B(TL{) the Banach space of trace-class operators on Wf with the trace-norm 
||A||i := tr^ f \A\, the Moyal product restricts to maps 

# : S kl (e, Ji(Hf)) x S k2 (e,B(Hf)) -> 5 fcl+fe2 (e, J^Ht)) 

and 

# : S k ^e, B(Ht)) x S k *(e, Jx{Us)) -> S k ^{e, MUf)) . 

The continuity of the product and of the asymptotic expansion holds as well with respect to the 
trace- norm || • ||i on Ji(Hf). 

One can compute the trace of a product of Weyl operators by integrating the trace of the 
point-wise product of the symbols: For scalar symbols b € S°(C) and a S W oa,1 (M. 2n ) = {/ G 
C°° (R 2n )\d a f e L 1 (M 2 ™) for all a G Ng™} it is known, e.g. |Ron84| . that a is trace class and that 

tr i 2 (E n)(aS) = 1 / dqdp a(q,p)b(q,p) . 



Now let A G W 00,1 (M 2n ,B(Hf)) and B e S°(B{H f )) such that tr Wf (AB) £ L^R 2 ™). Let (^i)jeN 
and (^fe)fc G N be orthonormal bases of L 2 (R") and respectively. Then au ■= {ifk, A(pi)u { € 
W oo,i( R 2n) and and ;= ^ uBipk ) Hi e 5-0( C ) for all fc, / g N. Thus 

j,k j,k,l 

= 77. 1 tz [ dqdp a k i(q,p)bi k (q,p) 
kt ( 27T£ > J 

= (2^ / d( l d P tT -Hf( A (<l,P) B (<l,P)) < °°, (9) 

where one should read this computation backwards to see that all expressions are well defined and 
that the equalities hold. 



2.3 The almost invariant subspace 

With an isolated eigenvalue eo there is associated a subspace WH of the state space that is almost 
invariant under the action of H. 

Proposition 1 (Super-adiabatic projection). Let Assumptions 1 and 2 hold. Then eo G 5°(M) 
and 7To G ( J7i (7^f )) . For e small enough, there exists an orthogonal projection TL £ G B(T~L) such 
that 

|| [iJ, IT] || = 0{e°°). (10) 

It is related to the band eo through 

||n e -7r|| = 0( £ °°), (11) 

where tt G S°(s, Ji{H t )) is a classical symbol with principal symbol itq taking values in the self- 
adjoint trace-class operators on Hf. H 6 is the spectral projection of tt for its spectrum near one, i.e. 
IF = Xji !](tt). The relevant diagonal block of the subprincipal symbol tti of it is given by 

TToTTlTTO = ^ 7r {7r , 7T }7T . (12) 

It is quite remarkable that for the following analysis we only need to know that such a projection 
W = 7r + 0(e) exists and that the diagonal block of its subprincipal symbol is given by (IT2"1) . 

With tt G S°(e, B(U{)) the statement of Proposition Q] was proven several times in the literature 
under slightly varying technical assumptions. The strategy of proof is due Nenciu and Sordoni 
[NS04] based on work of Hclffer and Sjostrand IHS90J . The exact statement of Proposition Q] 
follows, for example, from inspecting Theorem 3.2 and its proof in |Teu03j using the statements in 
Section l2~2l concerning Moyal expansion in S k (e, Ji{%f)) and, as a starting point, ttq G S°(Ji(Hf)). 
To see the latter, note that ttq takes values in the rank one projections. Thus by induction starting 
with 

dj-Ko = TT {djTTo)^o + ^o(9]^o)^o (13) 
all derivatives 0"ttq are of finite rank < N a and therefore trace class. Hence the fact that tt G 
S°(e,B(H { )) implies by 

sup ||3"7r (z)||l = sup tr Wf |9"tt (z)| < N a sup ||9^7r (z)|| < oo 

zGR 2 ™ zGR 2 ™ zGR 2 ™ 

for all a G N 2 ™ also tt G S Q (Ji{U t )). 
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Corollary 1. Under the hypotheses of Proposition\j]the diagonal Hamiltonian IP-HTP +n iZTI 
is self-adjoint on D(H), the projection Ii e almost commutes with the unitary time evolution operator 
in the sense that 



-iHs 



IT 



0(e c 



(14) 



iHs 



i(n e j?n [ +n e± ffn e± )i 



= 0(e°°\s\) 



(15) 



Moreover, for any A £ S {B{%f)) with k > 2n + 1 the operators Att and All e are trace class with 

tru(ATr) = 0(e- n \\A\\ L i) and tv n (All e ) = 0(e- n \\A\\ Ll ) . (16) 

Proof. The self-adjointness of the projected Hamiltonian on D(H) follows from the fact that [H, H e ] 
is a bounded operator. A standard Duhamel expansion yields (fl4l and (fTS"]) . For the last statement 
note that by © the estimate 

/ dqdp \tr Ul (A(q,p)w(q,p))\ < / dqdp |tr % (7r(g,p))| \\A(q,p)\\ <C\\A\\ L i < oo 

implies An £ J7i(H). Since IF is the spectral projection of 7r to its spectrum near one, there is a 
bounded continuous function g : M. — s- M. with \\g(Tt) — 1-h\\ = 0(e°°) such that IF = g(fc)n. This 
implies that also A H s is trace class and that 



tr H (irF) = tT H (Ag(n)7t) - tv n (An)(l + 0(e°°)) 



□ 



Note that (|15j) is the starting point for an approximation based on the concept of effective 
Hamiltonians. The idea is to map the range of IF unitarily to i 2 (R") using an eigenfunction 
ipo(q,p) in the range of iro(q,p). The image of the block IFflFF of H is then a pseudo-differential 
operator h e s, the effective Hamiltonian. Its scalar symbol h c s has an asymptotic expansion that 
can be, in principle, computed to any order. But it depends on a gauge, namely the choice of 
<Po(q,p)- For details on this approach we refer to HS901 [Teu03| . 



3 The classical Hamiltonian system 

For the semiclassical analysis within the subspace U e H we can now focus on the restricted Hamil- 
tonian H £ HH £ . The central object is the following scalar Hamiltonian function 



h(z) := e (z) +£,-z + e tr« f (H^z)^)) + sM(z), 



where 



M := |tr Wf ({ttoIFoKo}) . 
We will see that M can be written equivalently as 

M = |tr % ({tt q \H - e Q \ir }) = \ tr Wf ({tt q \H - eo|7r }vro) 
= -5 tr Wf (7r {7r , H - e }) , 



(17) 
(18) 

(19) 
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which is the expression found in |EW96j . who explain how to relate it also to [LF91J . The scalar 
Hamiltonian h has the property that its quantization approximates the action of the full operator 
H on the range of the projection IF up to terms of order e 2 , 

rF#iF = n e hn e + o{e 2 ). 

Proposition 2. The symbol h(z) :— h{z) — £ ■ z belongs to S (s, C) and h satisfies 

7T#/l#7T - 7T#iJ#7T = C(£ 2 ) in S° (g, B(U{ )) . (20) 

As a consequence, 



if/iIF - if tin 6 



0(e 2 ). (21) 



Note that, in general, there is no scalar symbol h such that (j2"0l) holds with an error 0(e 3 ). 
Neither does an analogous scalar symbol h exist in the case of a degenerate eigenvalue eo. Also note 
that in the introduction we did not distinguish between H and its principal symbol Hq. Indeed, 
if we replace Hq in all the following by H and ignore terms containing Hi, all the results hold 
equally. Of course then eo and tto need to be eigenvalues and spectral projections of H. In concrete 
applications, however, it is typically easier to split off the principal symbol. 

Proof. The claim that h £ S°(e,C) is an immediate consequence of the assumptions. Since the 
linear part £ • z appears also as a scalar in H, it suffices to check (f2"U| for h and H = Hq + eH\. 
According to (|6]) the asymptotic expansion of nffHffn starts with 

7r#£f#7r = eo7ro + e(7roiTi7ro + eo7ri7ro + eo7r 7ri 

- i7T {tfo,7ro} - ffo}7ro - 3N I #o \n }) + 0(e 2 ) 
= e 7To + e(tr« f (n Hi) ttq + e 7ri7ro + eo7r 7ri 

- %ir {H , M ~ |{tt , H }tt - |{7r |i/ |7r }) + 0{e 2 ) . 

On the other hand 

7r#/i#7r = e 7r + e(tr-H f (ttqHi) ttq + e 7ri7r + e 7r 7ri 

+ A/7T - 57r {eo,7r } - ^{710, e }7r - 5{vro|eo|7r }) + 0{e 2 ) . 

It remains to show that 

jMttq = -n {H - e ,ir } - {ir Ql H - e Q }ir - {ir Q \H - e Q \n Q } . (22) 
To see this note that §5§ and (IT5|) imply 

Mrt = ~tr nt ({7r |i?oko})7r ^ ±tr K( ({vr |i?o - e |7r }) tt 

= ^tr^ (7r {7ro|-ffo — eo|7r }7r ) 7r = ^7r {7r |77 — e |7r }7ro 
= %{tt \H - e |7r } . 



Moreover 



= dj ((H - eo)7r )7r = dj{H - e )7r + (H - e )djir iro 
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implies that 

{no, H - eojvro = -{tt \H q - e |7ro}vro = -{ir \H - e |7r } , 

which shows also (|19[) . By the same reasoning we have also 7To{-ffo — eo,7To} = — {ttq\Ho — eo|7To} 
and thus J221) follows. □ 



In order to obtain semiclassical approximations up to errors of order s 2 , we need to take into 
account that the restriction to the range of IP also induces a modified symplectic form u> £ on R 2 ™ 
given by 

( o E n \ ( n qq n qp \ 

u> £ := u +sSl =: ^ ^ Q )+£^ nPq QPP ) , (23) 

where the components of Q in the canonical basis are, with z = (q,p) and a, j3 = 1, ... , 2n, 

Ma/3 ■= -itr-H f (7ro[9 ZQ 7r ,9 Z|9 7r ]) . 

As is well known and will be shown also in Proposition [6j is the curvature 2-form of the Berry 
connection. 

The Liouville measure A £ associated with the symplectic form u> e has the expansion 

A e = (l + ietr« f ( 7 r {7r , 7 r }) + O(e 2 ))d (? 1 A---Adp". (24) 

We postpone a discussion of these objects to Section^ There we prove that w e is indeed a symplectic 
form and that its Liouville measure is given by (|24p . First, however, we present and prove our main 
results concerning semiclassical approximations. 

4 Semiclassical approximations 

We now explain how to approximate quantum mechanical expectation values 

tr-H (pa) 

of semiclassical observables a := a <S> 1% for the slow degrees of freedom in terms of the classical 
Hamiltonian system described in the preceding section. Of course this can only work for states p 
that belong to the range of n e . 

4.1 Expectation values for stationary states 

In this section we consider stationary states that are suitable functions of the Hamiltonian, i.e. 

P = f(H). 

Since we make use of the Helffer-Sjostrand formula for analyzing functions of self-adjoint operators, 
a natural space of admissible functions / is the space 

A := {/ 6 C°°(K) | 3/3 < : sup |(a;)'^' 3 f (n) (x)\ < oo for all n G N } . 

The following theorem shows that we can approximate the quantum mechanical expectation values 
by a classical phase space average up to errors of order e 2 if we use the correct classical Hamiltonian 
h and the correct Liouville measure A e . An important application of this result is the possibility to 
express quantum equilibrium distributions in terms of the classical system. 
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Theorem 1 (Equilibrium distributions). Let Assumptions 1 and 2 of Section \2.1\ hold. Then 
{or f € A and a £ S k (C) with k > 2n + 1 it holds that 

tr H (ll s f(H)a) = — L- (J dA e f(h(q,p)) a(q,p) + O (e 2 \\a\\ L ,] 

Note that the precise error estimate in terms of 1 1 a 1 1 ^ i is very useful when taking e.g. a thermo- 
dynamic limit, i.e. when looking at a sequence of observables a n that extend over larger and larger 
regions of phase space. This is done explicitly in Theorem [3l where we use Theorem [T] to compute 
the free energy per unit volume in the Hofstadter model. In that model it is also explicit, that 
only the phase space volume related to each band is modified, while the total phase space volume 
is unchanged as the total integrated curvature of all bands is zero. 

Proof. According to CorollaryHJ IF a is trace class with tr-^IFa) = C(£~ n ||a||Li) and thus we can 
estimate expressions of the form tr-H(RaH £ ) with R € by Ce~"||i?|| ||a|| L i with a constant 

C independent of e and R. Recall that \\[H , IF]| = 0(e°°) and, since h is scalar, \\[h, IF]|| = 
\\[h, tt]\\ + 0{e°°) = O(e). In the following computation we use these observations together with 
the statements of Lemma [T] and Lemma [2] from the appendix where indicated, 

tv n ^f(H)aU e ) # tr„(n £ /(#)II £ a) + 0(e°°\\a\\ L1 ) 

IP tr n (wf(U s HW)Wa) + 0{e°°\\a\\ L i) 



^ tr„(rF/(iFMF)rFaJ + 0(e 2 -"|H| L1 ) 

IP tr H (n e /Wn £ a) + 0{e 2 - n \\a\W) 
= ti n Uf(h)Tta) + 0(e 2 - n \\a\\ Ll ) . 



Next note that for scalar symbols the functional calculus for pseudo-differential operators implies 
that \\f(h) — f(h)\\ = 0(e 2 ), c.f. e.g. Chapter 8 in [DS99] . Since this is an important and nontrivial 
step, we sketch a proof of this statement in Lemma |3] in the appendix. Hence with fh :=/o/i£ 
S° (e, E) we have 

tv n (ll E f(H)a) = tr« (nf(h)™) + 0{e 2 - n \\a\\ L i) (25) 
@ 1 



(2?re) n 
1 

(271-e)" 



dqdp tr n{ ((n#f h #ir)(q,p)a(q,p)) + 0(e 2 - n \\a\\ Ll ) 
dqdpa(q,p)tT n[ ((TT#f h #ir)(q,p)) + 0(s 2 - n \\a\\ L r) , 



where 7r#/ h #7r e S°(e, Ji(U f )). 

It remains to compute tr% { (7r#//j#7r) up to 0(e 2 ) in S°(e, Ji(Hf)). Using ((6]) and the fact that 
fh is scalar, we have that the expansion of F :— Tr#fh#Tr starts with Fq — irofha and 

Fl = TTofhl^O + fhOKiir + fhO^O^l - h^o{fhO, TTo} - -H^O, fho}^0 - hfho{^Q, 7To} ■ 
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Taking the trace we get 

tr Wf (7r#/ h #7r) = ix Ht (ttoAo + eF x ) + 0(e 2 ) 

= /fco + e(tr« f (7r -Fi7r )+tr Wf (t^F^)) + 0{e 2 ) 
= fho + efhi +e2f h0 tr H{ (ttottitto) + 0(e 2 ) , 

where we used that for scalar symbols 

{fh,Ko} = 7ro{/h,,7To}7r +7T Q {/h,7ro}7To 

has vanishing trace and that also tr^ f . {ttq, ttq} = 0, see With the expression (|12p for i\\ the 
final result is 

tr« f (7r#/ fc #7r) = f h (1 + ietr Kf M^o, 7r })) + £>(e 2 ) . 
Inserting this into (l2"5l and comparing with (|24l) completes the proof. □ 

4.2 Egorov theorem and transport of Wigner functions 

For non-stationary states p(t) — e~ lH i pe lH i starting in the range of IF, i.e. p = IF pip, (TT41 
implies that 

U £ p(t)U £ = Il e e- i£[ e pe i6 * n £ = e" 1 ^ WpSFd*' + 0(\t\e°°) 
= p(*) + 0(|t| e ~). 

For semiclassical observables a = a ® l« f we thus obtain 

tr«(p(i)a) = tr^p IT e 1 ^ ae^f* IF) + Q(|t|g°°) . 

= :A(i) 

The Egorov theorem shows that on the range of IF one can approximate the quantum mechanical 
evolution A(t) by evolving the symbol a along the classical flow </>* up to errors of order e 2 . Its 
proof for scalar Hamiltonians follows from the relation ((5J) between commutators of operators and 
the Poisson bracket for the symbols. As the following proposition shows, a similar statement holds 
in our case, but with a modified Poisson bracket. 

Proposition 3. Let Assumptions 1 and 2 of Section \2.1\ hold. Then for e small enough, the 
Hamiltonian vector field Xf t = — (u> e ) al3 dph is a classical symbol in S°(e,M. 2n ) . There is a r G N 
such that for any a G S (e, C) it holds that 

n s (i[H,a] - Op w (AWa))n e | = O(e 2 \\a\\ ,r) ■ (26) 

Here Xh ■ Vet = {h,a} Ue is just the Poisson bracket with respect to the modified symplectic form. 
For a(z) — r] ■ z with r\ G R 2n formula i26\) holds with an error of order 0(£ 2 |?7|). 

Proof. The statement about X h follows immediately from n a/3 G S°(R) and V/i G S°(s,M 2n ). Let 
a G S°(e, C), then with dTOj) we have 

in £ [i?,a]n £ = i[n e OT E ,n e an E ] + o(e°°||d||). 
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In order to replace H by h, observe that 

ifffj - h)IT = iOp w (7r#i7#7r - 7T#/i#7r) + 0(e°°) 

and that by M 



ha := s- 2 {it#H#it - 7r#h#7r) G S°(e,B(%)) 
satisfies 7r#ft,2#7r = + 0(£°°). Hence its principal symbol satisfies 

(Mo = ^0 (h 2 )o TO 

and we obtain for sufficiently large r € N with r > 2n + 1 that 



n £ (ff-/i)n e ,irair 



1£ 



/lo, 7T a 7T 



•o(e°°l|a| 



= ieOp w ([(/i 2 )o, vroa^o]) + O^IMIo.r) - O(e 2 ||a||o,r) 
The remaining term can be rearranged in the following way, 



ir/inMPair 



= ±W[h,a\W + iff ([h,H eX ] [a,U e± ] - [a,U s± ] [h,U £± ]^j IT 

= irr [h, a]n e - k it U [h, it] , ± [a, it] 1 it . 

Since ft. and a are scalar, © implies 

i[M] - Op w ({/ l ,a}) + 0( £ 2 ||a||o, r ) 
i[/i,n s ] = Op w ({ft,7r }) + O(e) 
i[a,lT] = Op w ({a,^o}) + O(e||a||o, r ). 

Here and in the following steps a priori different rs might be necessary. We agree to denote by r 
the largest one that works in all estimates. Summing up the estimates we got up to now we have 



^W[H,a]W = IT Op w ({/i,a}-i£[{/i,7r },{a,7r }]) IT + 0(e 2 



Wo,- 



The subprincipal symbol on the right hand side can be replaced by a scalar one at leading order, 
because of the outside projections, 

IT Op w {[{h, 7T }, {a, 7T }]) IT = n £ Op w ( 7 r [{/ 1 ,^ },{a,^o}]^o) IP + O(e||a|| ,r) 
= n £ Op w (tr^(7ro[{h,7ro},{a,7ro}])7ro) IT + O(e||o|| 0l r) 
= n e Op w (tr Wf K[{/ l ,7r },{a,7r }]))n e + O( £ ||a|| , I .) 

and thus 



in E [if,a]n e = it op w ({M}-ietr« f (7ro[{Mo},{a,7ro}])) n £ + o( £ 2 



mo,. 
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Comparing this with 

(Va) • X h = (d q a, d p a) 



-e WP E n + eQ pq 



d q h 
d p h 



+ o(e 2 lkllo,i) 



= {h, a} - ietr % (ir [{h, tt }, {a, vr }]) + O(e 2 ||a|| ,i) 

completes the proof for a 6 S°(e,C). Now notice that the result implies that on the level of classical 
symbols we have 

^#(H#a-a#H)#n = 7r#(X h -Va)#7r + 0(e 2 ). (27) 

However, as an algebraic relation it holds for a in any of the symbols classes S k (M), in particular also 
for a(z) = Tj-z. But then ±(H#a-a#H) = {H,a} e S°(e,£(U{)) and (X h -Va) = X h -Tj £ 5°(e,M) 
imply that also the remainder term is C(e 2 ) in S°(e, C(Hf)). 

Note that checking relation (j2"T)l directly by expanding both sides in powers of e is an alternative 
way to prove the proposition. This was done in [Still] and leads to much more painful computations 
than the argument given here. □ 

Theorem 2 (Egorov theorem). Let Assumptions 1 and 2 of Section \2.1\ hold. Then the Hamil- 
tonian flow <j)\ of (h, ui e ) exists globally and for any a s 

S°{e,C) H holds that 



IT (A{t) - ao ft 
uniformly on bounded time intervals, where 

A(t) := c [ "i 



)n e = 0(e 2 ) 



,-iff4 



Proof. Since the Hamiltonian vector field Xh is smooth with all derivatives uniformly bounded, the 
corresponding flow (^*) tg R exists globally. Each map 0* : R. 2n — > R 2 ™ is a C^-diffeomorphism with 
all derivatives bounded. Hence t n- ait) := a o 0* g S (e, C) is smooth and uniformly bounded on 
bounded intervals in time. In order to apply Proposition [3] we use the standard Duhamel argument 
and (fill) to show that 



IT (A(t) - o(t)) IT = n £ J ds-^ (e iA ia(t 



-iHi 



n- 1 



dse 1 



n e 



ds e 



H,a(t - s) 



^n e a(t- s )n E e 



0(e°°t 2 sup |[a(s) 
«e[o,t] 



Hence the claim follows from ([26]) and the fact that a(t) := a o 0* satisfies 4;0,(i) = Xh ■ Va(<). □ 

As a corollary we get the following result on the transport of Wigner functions, which is dual 
to the Egorov theorem. For i\> € L 2 (R",Hf) the Wigner transform of ip is defined by 



(2tt)" 



dx&»*ll> (q - fx) ® ^* ( 9 + fx) 
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For := \ip)(i/j\ and A G S°(B(H f )) it holds that 

tr H (p*A) = f dgdptr« f (W^(q,p)A(q, p)) . 

JR 2n 

Let u e be the density of X £ with respect to Lebesgue measure, i.e. A e = v e dqdp. Then for tpo G H £ TL 
and i/j(t) := c~ lH zip we find for any a G S°(C) that 

tr«(^ W a) = tr w (^M(t)) = tx n {p^W A{t)W) 

= tr„(^IF ^)n e ) + 0(e 2 ) = tv n (p^^t)) + 0(e 2 ) 

= [ dqdpty n{ (W*(g,p)oo^( ?l p)) + 0(e 2 ) 

= / dA £ Kfep))" 1 tr„ f (W^°( ?> p)) a o #( ?>p ) + 0(e 2 ) 
Jm 2 ™ <• v ' 

-■wt°(q,p) 

= / dA e <° o <fc\q,p) a(q,p) + 0(e 2 ) , 

where we used that the Liouville measure A e is invariant under the Hamiltonian flow <j> e . Thus we 
define the reduced or scalar Wigner transform wf associated to the projection tt by 

wt := {v s {q,p))-hv Ul W^ = (1 - ietr % (7r {^ , M)) tr W( W^ + G(e 2 ) . 
Corollary 2 (Transport of the Wigner function). For ipo G H £ H it holds that 





w m 

WO 


= W^ ° o <; 


i>7* + 0{e 2 ) 


in the sense that with ip(t) := e lH 


^0 






(V>(t), &1>(t)) H 


-/ 

«/]Ri" 




(q,p)a(q,p) 






t dA £ «° 


o ( f>7 t )(q,p)a(q,p) + 0(e 2 ) 



for all at 5°(C). 

5 The geometry of the classical phase space 

In this section we show that u e is a symplectic form, compute its Liouville measure A e and clarify 
its relation to the geometry of the line bundle defined by the family of projections TT (q,p) over R 2 ™. 

Proposition 4. For e small enough u s defined in \23\) is a symplectic form. 

Proof. The matrix Q is by definition skew-symmetric and bounded, hence detw £ ^ for e small 
enough. Thus uj £ is a non-degenerate 2-form. To see that ui e is closed, note that 

i<9 27 f2 Q/3 = tr-Hj (dg y Tro[dz a iro,d Xe iro]) 

+ tr« f (TT [d Zy d Za 7ro,9 Sj8 7ro]) + tr Hf (ir [d Za Tr Q ,d Zy d Zf3 ir }) 
= tr K( (7ro[^9 ZQ 7r ,5 Z|9 7r ]) - tr« f (tt [d^ 5^ tt , 9 ZQ 7r ]) 
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implies 9 ZT f2 Q ^ + d Za fl^ 4- d Zft Q-y a — and thus 

d (n a pdz a A dz' 3 ) = d Zy n aP dz~* A dz a A dz 13 = . 

□ 

Proposition 5. TTie Liouville measure of the symplectic form uj £ is given by 

X E = (l + ietr % (7r { 7 r ,7r }) + O(e 2 ))d (? 1 A-.-Adp". (28) 

infn— 1) 

Proof. The Liouville measure is denned by A e = — — — , u>£ n . With this sign convention Wq" 

yields the canonically oriented volume form on R 2 ". To compute the first order corrections, it 
suffices to take all terms into account which contain exactly one term of the form £l a p dz a A dz@ . 
Due to the antisymmetry of the wedge product, only the off-diagonal terms f2^? dq 3 A dp 3 and 
SVjj dp 3 A dq 3 contribute and we get 

A e = f 1 + | J2 - Aft) j dg 1 A • • • A dg" A dp 1 A • • • A dp" + G(e 2 ) 
= (1 + ietr Hf (n {n , tt })) dg 1 A ■ • • A dp" + C(e 2 ) . 

□ 

We already mentioned and it is well known that ft is, up to a factor of i, the curvature of the 
Berry connection. For the sake of completeness we explain what exactly this means. Originally the 
Berry phase was introduced for the adiabatic limit of time-dependent Hamiltonians by Berry |Ber84] . 
Shortly thereafter Simon [Sim83j realized that this phase could be rewritten as the holonomy of the 
curvature of a certain line bundle. Let us briefly explain this idea using our notation. The trivial 
Hilbert bundle E := M. 2n x Hi A M 2 ", where P is the projection onto the first component, 
is endowed with a canonical flat connection: for -0 G T(i?) and X G L(TR 2 ") let (V ' x^){z) = 
X 3 d z .ip(z). Within this picture we can regard the symbol Hq : M. 2n — > B(Hf) as a section in the 
cndomorphism bundle of E, i.e., Ho G L(End(i?)) which acts on sections ip G Associated to an 

isolated eigenvalue eo : R 2 " — > M of H is the spectral projection n : R 2 " — > B(7i{) which again we 
regard as a section ttq G L(End(_E)). This projection defines a sub-vectorbundle L := {(z, ip) | ip(z) G 
7To(2:)'Hf} of E. Now the trivial connection on E induces a connection on L by projection, i.e. 
V x 4> '■= ttoV x<t> for <f> G r(L) C T(E). This connection is called the Berry connection. 

Proposition 6. The curvature form R' of the Berry connection V is ^R'^dz 1 Adz 3 , where 

R'ij = 7r [9 2i 7ro, <9 Zj .7r ] 

and for a line bundle 

R'ij = tr Wf {■Ka[d Zi 'K Q ,d Zj 'K ]) . 
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Proof. Let X,Y e r(TIR 2 ") and <fi e T{L), then by definition 

R'{X,Y)ct> := (VxVy-VyV'x-V^y^ 

= 7T Vx (7T O Vy0) - 7T Vy (7T O Vx0) ~ 7T O V[X,Y]0 
= (TTQ^ (7109,(7100)) - TTofli (^^(TTQ^))) 

= X'F- 7 (7ro[9i,7ro]9j(7r o 0) -7ro[9,-,7ro]9i(7r ^)) 

= ^^TTopiTTo.fliTro]^. 

Here we used 7To</> = <f> for G T(L) in the third equality, commutativity of partial derivatives in 
the fourth equality and (|13p in the fifth equality. If rank7io = 1, then 

7r [9i7ro, djTT ](j) = 7r tr Wf (7r [9i7r , dj^o])<j> = tr« ( (7ro[9i7To,9,-7ro])^ • 

□ 

Since R 2n is contractible, there exists a global section (ipi, . . . , ipd) of the frame bundle associated 
with L, i.e. tp±, . . . ,ipd £ T(L) and (ipi(z), . . . , tpd(z)) is an orthonormal basis of TTo(z)Hf for each 
z € M 2 ". With the physicists' bracket notation the projectors can be written more explicitly in the 
form 7To(z) = \ife(z)) (ipi(z)\ and R'^ can be expressed in terms of this frame. For the case d — 1 of 
a line bundle one finds 

R'ij = {diV{z),djLp{z)) - (d j ip(z),d t ip(z)) . 
The connection one-form A of the Berry connection with respect to this trivialization is given by 

A = Ajdz j = -i ((p(z),djip(z)) dz j , (29) 

i.e. it acts on a local section cj> = ipip as 

V dj <f>=:V d .(^) = ((d j +iA j )i>)<p. 

The curvature two-form Q is just 

il = dA, i.e. ftij — diAj—djAi. 

We remark that the Berry connection is usually discussed in the form of the connection coefficient 
(|2^|) with respect to a local or global trivialization. Note that for our discussion of the semiclassical 
approximations we used at no point the existence of a trivialization of the bundle L. This is why the 
same reasoning can be used also in the case of Bloch electrons in magnetic fields. There the phase 
space is the cotangent bundle of a flat torus and the corresponding bundle L is not trivializable, 
c.f. Section [7] 

6 Discussion of the literature 

Semiclassical approximations for Hamiltonians with operator valued symbols have been discussed 
many times in the literature and almost all objects appearing in our analysis appeared in some form 
before. Hence we give a short roadmap to the literature that we know. 
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Our first reference for the corrected symplectic form is a series of papers by Iida and Kuratsuji 
KI85] IIK871 IKI88] . They used path integrals to study adiabatic decoupling of a slow and a fast 
system. In the path integral formalism they derive an effective Lagrangian for the slow system that 
incorporates Berry phase effects. After taking the semiclassical limit they are able to write down 
their version of the semiclassical equations of motion. However they missed the M term (|18p and 
thus their result must be considered incomplete. 

In the context of WKB approximations the work of Bernstein |Ber75j is the earliest reference 
which contains an additional phase term besides the Berry phase. Kaufman, Ye and Hui [KYH87J 
used a variational approach to re-derive those results and expressed them in a form using Poisson 
brackets. The comparison of these results to the more modern ones is not straightforward. An 
attempt to reconcile them has been made by Fukui [Fuk92 . 

A thorough discussion of the older literature and the first appearance of the general M term 
and the modified symplectic form lo £ is contained in the work of Littlejohn and Flynn |LF91j . They 
use the Weyl calculus to diagonalize operators with matrix-valued symbols. To do this, they must 
assume that the spectrum consists only of globally isolated eigenvalues and that the linebundle L 
defined in the previous section is trivializable. They found that - unlike in the scalar case - the 
equation for the amplitude in the WKB approximation contains additional phase terms. One of 
them is the Berry phase and the other the M term, which they call "no-name term" . The Berry term 
is gauge dependent, i.e. it depends on the choice of a global tirvialization. This gauge dependence 
seemed to be unsatisfactory from their point of view and they searched for a method to express 
their results in a gauge-independent way, which turned out to be closely related to the right choice 
of coordinates. A few lines of calculation show that the phase space coordinates z'j := Zj — eAj(z) 
are Darboux or canonical coordinates. Here Aj are the components of the Berry connection (|29|) . 
This simply means that in these coordinates the symplectic form u e is given by the canonical form 
dq' A dp'. Littlejohn and Flynn on the other hand started out by canonical coordinates and the 
analysis of gauge invariance led them to consider non-canonical coordinates, which are associated to 
the symplectic form u e . However it seems to us, as if they failed to realize that their non-canonical 
coordinates are in fact position and momentum of the slow degrees of freedom. At least they neither 
proved nor claimed the statements of our Theorems [1] and [2j 

The algebraic part of the construction of the Moyal-projection ir for £-pseudo-differential op- 
erators was first done by Helffer and Sjostrand HS90 . Independently and motivated by |LF91j . 
Emmrich and Weinstein |EW96| gave a derivation of the transport equation in WKB approxima- 
tions which does not rely on eigenfunctions and is intrinsically gauge-independent. In their work 
they give an alternativ construction of the Moyal-projection ir and find the gauge independent 
classical Hamiltonian h with the M term in the form (|19| . For them the modified symplectic form 
and the modified flow <j)\ play no role. 

Completely independently, a formal derivation of the modified semiclassical equations of motion 
(fTj) in the special case of Bloch electrons using propagation of wave packets was given by Chang 
and Niu in [CN96j . 

Based on the approach of [HS90], the general construction of the super-adiabatic projection 
n e was first done by Nenciu and Sordoni NS04J. Using an adaption of this construction and 
in addition global trivializations of L, in PST03] the first rigorous derivation of the modified 
semiclassical equations for Bloch electrons without strong magnetic fields was given. The fact that 
such a trivialization is not available in magnetic Bloch bands was the motivation for the present 
work. 

There have been many more works considering the semiclassical limit of particles with spin, 
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which corresponds to Hf — C N . Among them we mention |BGK01[ IBG04] , where in |BG04j an 
Egorov theorem to any order in e is proved, but not based on a modified classical flow. 

To our best knowledge the statement of Theorem Q] on stationary expectations was used only 
recently by Gat and Avron |GA03] and later in [XSN05] . As far as we understand, there the state- 
ment is taken for granted, but not systematically derived or proved. The statement of Theorem [5] 
was proved in a special case in |PST03) . but otherwise neither claimed nor proved. 

Finally we would like to mention that the physics literature on the semi-classical limit for 
operator valued symbols is vast and it seems that many observations are reproduced independently 
several times. For example in [GBMQ7] ideas appearing in |LF91j are reproduced without citation. 
Since all these newer works we know of are at most formal, make no statements similar to ours and 
differ from our strategy of proof, we do not list them here. 

7 Application: Magnetic Bloch bands in the Hofstadter model 

Consider a gas of noninteracting fermionic particles on the lattice Z 2 subject to a constant magnetic 
field B = ( _°b % ) with Bgl. The single particle Hamiltonian is the discrete magnetic Laplacian 

H B = £ T* 

|n|=l 

acting as a bounded self-adjoint operator on £ 2 (Z 2 ). Here (T B ip)j = e^ n ' B: 'ijjj-n for n G Z 2 are 
magnetic translations. This is a simple model for particles in a periodic potential in the tight binding 
approximation. For B = the operator H B=0 is invariant under lattice translations and thus 
diagonalized by the Fourier transformation T : t 2 (I?) — > L 2 (R 2 /(2ttZ 2 )), where FH B=0 F* = E 
becomes the multiplication operator with the function E(k) = 2(cos(fci) + cos(A;2))- 

For B ^ 0, the operator H B is invariant under the dual magnetic translations (T B ip)j = 
e~^ n ^tpj_ n , which, in general, do not form a unitary group representation, since T^^T^^ = 

c lB TR iyT^ Q y However, for magnetic fields B = Bq = 2n | that are rational multiples of 2n, the 
magnetic field and the lattice are commensurable in the sense that the magnetic flux per unit cell 
is a rational multiple of the flux-quantum $o = = 27r. As a consequence one can extend the 
dual magnetic translations to a unitary representation of the subgroup T q := {7 G Z 2 | 71 6 qZ} of 
Z 2 by defining 

h^Bq / rpBo \ f rpBo \ 

J 7 — V J (l,0)j V C ' 1 )/ 

and thus H B ° can be fibered with respect to this group representation: Let M q :— [0, 2Tr/q) x [0, 2ir) 
be the reduced Brillouin zone and define the magnetic Bloch-Floquet transformation as 

U Bo :£ 2 (Z 2 ) ^ L 2 (M q ;C q ), 

{i>j)j^ -> (U B °iP) m (k) := e- ik - m e ifc ^(Tf^) m form = 0,...,g-l. 

Then U B ° is unitary and because of {U B °il>) m {k + 7*) = e~ i ' 1 "'- m {U Bo ijj) m (k) for any 7* in the dual 
lattice T* =:= {(71,72) € K 2 | 7i G yZ, 73 G 2-7rZ} the functions in its range can be naturally 
extended to r-equivariant functions in L 2 oc (IR 2 ; C 9 ). Here a function / : IR 2 — 5- C 9 is said to be 
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r-equivariant if it satisfies 

/(fc + 7*) = r( 7 *)/(fc) :=diag(l,e-^. 



<i-^i)f(k) for all 7* ef! 



for the unitary representation r of T* . 

The image U B °H Bo U B °* of the Hamiltonian H B ° acts on "H = L 2 (M q ; C q ) as the matrix- valued 
multiplication operator 



/ 2cos(fc2) 



H (k) 



2cos(fc 2 + Bo) 







2cos(/fc 2 + 2B ) 










V 



ifcl 



e ifcl 2cos(fe 2 + (?-l)S ) / 



which is again r-equivariant, i.e. Ho(k+j*) = r(7*)ifo(fc)'?"( — 7*) and thus r*-periodic up to unitary 
equivalence. It is well known that for q odd, Ho(k) has g distinct eigenvalues e^-'(fc) < . . . < e^(fc) 
such that e^'(M q ) (~]e^(M q ) = for j ^ i. By the r-equivariance and analyticity of Ho(k), e^\k) 
are r*-periodic real analytic functions. For q even, all bands but the middle ones are isolated, 
where e^ 2 )(A;*) = e( ? / 2+1 '(£;*) for one point £;* G M g . Thus except for the middle bands the same 
conclusions hold. The spectrum of H B ° thus consists of q resp. q — 1 disjoint closed intervals and if 
plotted as a function of B$ gives rise to the famous Hofstadter butterfly with its fractal structure. 

We now show that small (also irrational) perturbations of B give rise to a semiclassical problem 
for a Hamiltonian with operator valued symbol and that we can use the semiclassical approximations 
developed in the previous chapters in order to understand H B for a magnetic field B = B + eb in 
terms of a classical system based on the geometry of Hg(k). To see this note that for B = Bo+eb and 
in the presence of an electric potential V £ (n) = V(en) with V{x) = Vb(x) + £ ■ x for Vb G C^'(M 2 ), 
the Bloch-Floquet transformation with respect to Bq yields the Hamiltonian 



H £ := u Bo H B U Bo * = Hoih + \bied r h2 , k 2 



acting on L 2 (M q ; C q ). Here 9 r . denotes the derivatives with r-equivariant boundary conditions and 
the matrix entries, e.g. 2cos(fc2 — ^bied^), are defined by the functional calculus for self-adjoint 
operators. However, the operator H e can also be understood as the e-Weyl-quantization of the 
matrix-valued function 



H(k,r) := H (k! + \br 2 ,k 2 



V{r) 



i.e. H e = H(k,ieVJ,). In order to apply the Weyl-calculus from M 2 here, one identifies L 2 (M q ,C q ) 
with the space L 2 := {/ 6 i 2 oc (M 2 ) | / is r-equivariant} with the norm \\f\\ 2 = -^p- J M \f{k)\ 2 dk 
and restricts to T-equivariant symbols, see [PST03] or Appendix B in |Teu03| for details. With this 
modification all the results of this paper hold with the identical proofs for the Hamiltonian with 
symbol H(k, r) and the underlying classical system on the phase space T 2 x R 2 , which we describe 
next. 

Let e^'(k) be an isolated band of Ho(k) with spectral projection 7r^(fc). Then e^\k,r) := 
(k+hbr)+V (r) is an isolated band for H{k, r) with spectral projection fr^ (fc, r) = 7r^ (fc+ibr). 
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A simple computation shows that the associated classical system is 

ft«) (k, r) = e (j) (Jfe + |br) + V(r) + sbM ij) (k + |br) , 

where 

= Im trr? 

is the effective magnetic moment of the particles in the jth band. With 

ft« := -itroi^^Tr^^Tr^]) 
being the curvature of the Berry connection on the jth eigenspace bundle of HqQc) one finds 

amk (k )_( o + ibr) \ y)r i 

" (fc ' r) - ^ _n(j)(fc + i b r) J-** ( fc +2 br )> 

q(j>* = -I b rj (j) = ibQ^E 2 and n<» r = -±bfi (i) b = i&Q^b. The corresponding equations of 
motion (JTJ) become much simpler if one changes coordinates to the kinetic momentum k :— k + ^br, 
leading to 



-b E 2 
-E 2 enW(«) 



h {j) {r,n) = e (j) ( K ) + V(r) +£&X (j) (k) and = 
with Liouville measure AXi = (1 + e&fi^ (k)) drdK. The Hamiltonian equations of motion are thus 



(l + eMlCj)(«)) ^ -£ 2 b ) { d K h( K ,r) ) 

1 / d K ( e W)(«) + e6A<W)(«;)) - en^(/c) d r y(r) \ 

(l + e6nW)(«)) ^ -5 r F(r)+b9 K (e«(/c) + J' 

or, equivalently, by 

r = d K (e {j) {K) + ebM {:j) {n)^ - eVL {j) {n) k, k= -d r V(r) + br . 

Finally note that by r-equivariance of Ho(k) the functions My* and QP> are r*-periodic. 
However, due to the symmetry of the original problem they are even periodic with respect to ^Z 2 . 
Therefore the classical system can be defined on the phase space T q x R 2 where T q is the torus 
[(W<Z) 2 - 

We now compute the free energy per unit area of a gas of noninteracting Hofstadter particles at 
temperature T — 1//3 and chemical potential fi, also called the pressure. To this end let \ n : R 2 — > 
[0, 1] be a sequence of "smooth characteristic functions" on A„ := [— n, n] 2 , i.e. Xn £ 5°(R) uniformly 
with suppx„ G A„ and XnU„_i = 1- With a slight abuse of notation we set A„ := HxnlU 1 - 



= 2tt2 

9 

and let B = B a + eb and V(r) = 0. Then for < /3 < oc, /i £ R and e > small enough, the limit 



Theorem 3 (The thermodynamic pressure in the Hofstadter model) . Let Bq = 2tt | with q odd 



p(B,0,n) =j9~ 1 f Um 4^ tr ^(z 2 ) (xn(ex)ln(l + e-^ HB -^)) 
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exists and is approximated by the classical expression 

p(B, 0, (j,) = ^ 1 { 1 + £bnU) ^) ln i 1 + e-^M-ri) dn + 0(e 2 ) . 
Note that j^|A„| is the area of the region in space to which X n (sx) restricts. 

Proof. The existence of the limit follows immediately from invariance of the Hamiltonian under 
magnetic translations. By unitarity of the Bloch-Floquet transformation and Theorem Q] we find 

tr, 2(z2) (xn{ex) ln (l + e -« ffS -">)) = tr w (*„(kV£) ln (l + e^**""))) 
= £> w (nW) Xn (kV3E)In(l + e-^-"))) 

= / (l + £&O«( K )) Xn ( r )l n (i + e -^ W) W-^)dKdr + 0(||x n |UO 

^ {2-Key J MqXM 2 V /V / 



= £t|H>/ (l + e6nW)(«))ln(l + e-« hW3 W-/*))d« + 0(|A„|). 

Now we can restrict the integration to T g to get a factor of g, divide by the area £ _2 |A„| and take 
the limit to prove the second statement. □ 

The thermodynamic pressure is the starting point for computing several physically relevant quan- 
tities like the magnetization M(B,/3,fj,) = 8bp{B, (3, /x) and the density p(B,f3,n) = d fJi p(B, f3, fi) 
and derivatives thereof. 

Corollary 3 (The magnetization of the Hofstadter model at rational B). Let Bq — 2ir | with q 
odd and V(r) = 0. Then for < j3 < 00, fj, 6 M it holds that 

M(B , (}, M ) = ^ £ jf (ffi, M (h^)M^ + i In (l + e-W 0) -">) Q«) d. , 
where fp^(E) := (1 + e^ E ~^) _:l is i/ie Fermi-Dirac distribution. 

Proof. The formula is obtained by taking the derivative of the semi-classical expression for p from 
Thcorcm[3l Noting that 8b — b~ 1 8 £ shows that the error term does not contribute to the derivative, 
since it is 0(e 2 ). □ 

This formula for the magnetization in solids is rather recent even on a heuristic level, see 
GA03, XSN05 . It was proved by a completely different approach as a special case of a more 
general formula in [SBT12 . Note however, that at least on a formal level the present derivation 
using the semiclassical model is much simpler. 

As a final application of the Egorov theorem or, more precisely, of Proposition^ we compute the 
current density in the Hofstadter model at zero temperature when a constant electric field £ £ M 2 
is applied and when fj, lies in a gap between e^™-* and e^ m+1 \ The corresponding state is given by 
IF = jyjLi ^ ano - localization to the region £ _1 A„ by 11^ := H E Xn(£x)Il e . 
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Theorem 4 (Current density in the Hofstadter model). Let Bo = 2ir ~ with q odd and let B = 

B + eb and V(r) = £-r. Let y, £ R be such that e (" l )(fc) < fi < e { - m+v >{k) for one m 6 {l,...,q} 
and all k £ T q . Then with W = Y^jLi ^ an d Hn := ^ £ Xn(£x)H £ it holds that 

j(B,£):= lim - T ±- tt P[P) (W n J £ ) = £ ± jfa £ / n^(«)d« + 0(e) , 

where J £ := ^[H B ,x] is the current operator. 

As j- L n^(Ac)dK £ Z is the Chern number of the jth magnetic sub-band at magnetic field B Q , 
the Hall current is quantized and independent of B in a neighborhood of Bo at leading order. Since 
this result is basically well known and proven in greater generality by other means, e.g. BES94 , 
we only give the simple computation based on the semiclassical model, but skip the mathematical 
details. However, the following can be easily made rigorous by using the results of the previous 
sections. On the other hand we believe that the strength of the semiclassical model lies in the 
fact that its application is very simple and leads in a straight-forward way to the correct physical 
expressions. 

According to J2BJ) we have that J = ±(H#r - r#H) £ S°(R) satisfies 

noo j t jiw = no - ) x Kri n^ + o(e 2 ) 

and thus 

m m 
$>« (Xnll« J,n«) = (Xnri j) *fc,r,nW) + 0(e 2 \A n \) 

j=i i=i 

— f ( \ 

= (2^fE J d ^ ) (1+ ^)( K) ) ( d « ( e(J)(K) + - E " (i) ^ ■ 

3_1 T,xl 2 

bmce dividing by -r|A n | and taking the thermodynamic limit yields 

m „ 

= j^p^/ d«(^(e«)(«) + e hMC7)(«))_ en 0-)(«)£) 

Here we used that the leading term is the integral of a gradient over a manifold without boundary 
and thus zero by Stokes' theorem. This is the well known text book statement that filled bands can- 
not contribute to the current. The standard semiclassical model without the first order corrections 
would imply exactly this. 

A Technical lemmas 

In the following proofs we apply the Helffer-Sjostrand formula for functions in A as developed in 
[Dav95j . 
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Lemma 1. Let A be a self-adjoint operator on H with domain D and P £ B{'H) an orthogonal 
projection such that PD C D and 

\\[A,P]\\<d. 

Then for any f £ A there is a constant C < oo such that 

\\[f(A),P}\\ < Cd (30) 

and 

\\Pf(A)P-f(PAP)\\ < Cd 2 . (31) 

Proof. The bound ([50)1 is an immediate consequence of the Helffer-Sjostrand formula. For (I5T|) first 
note that PAP is self-adjoint on U with domain PD®P^'H. Hence f(PAP) £ B(H) is well-defined 
by the functional calculus and it holds that f(PAP) = Pf(PAP)P. We next show that 

||P ((A z)- 1 - (PAP z)- 1 ) P\\ < . (32) 

This follows from taking norms in the identity 

P ((A - z)- 1 - (PAP - z)- 1 ) P = P(A - z)- x P - P(PAP - z)- x P 
= P(A - z)- l P ((PAP -z)-(A- zj) P(PAP - zY l P 

+ P(A - z)- 1 P(A - z)P(PAP - zY l P - P(PAP - z)- 1 P 
= P(A- z)- 1 [P,A}P(PAP - z)- 1 P 
= P(A - - P)[P,A}P(PAP - z)- 1 P 

= P[P,(A~ z)- 1 }[P,A}P(PAP - z)~ 1 P 
= P(A - z)-\A, P](A - z y 1 [P,A]P(PAP - zY 1 ? . 

Now let / be an almost analytic extension of / with \d s f(z)\ < c|Im(z)| 3 . By the Helffer-Sjostrand 
formula and (1321) 



\\Pf(A)P-Pf(PAP)P\\ = 

= \ [ d,f(z)P((A-z)- 1 -(PAP-z)- 1 )Pdxdy 
Jc 

< \ I \d- z f(z)\\\P((A-z)- 1 -(PAP^z)- 1 )P\\dxdy < C f d 2 . 
Jc 

□ 

Lemma 2. Let A,B be self-adjoint operators on % with domain D such that A — B extends to a 
bounded operator with \\A — B\\ < d. Then for any f £ A there is a constant C < oo such that 

\\f(A)-f(B)\\ < Cd. (33) 

Proof. Clearly 

||(,4 - z)- 1 — (B — z)- 1 !! = \\(A- z)-\B - A)(B - z)- l \\ < ^?L_ . (34) 
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Let / be an almost analytic extension of / with \dgf(z)\ < c|Im(z)| 2 . By The Helffcr-Sjostrand 
formula and (l34l) 



\\f(A)-f(B)\\ = 



- / d- z f(z) ((A - z)- 1 -{B- z)- 1 ) Axdy 



< - I \d- z f{z)\\\(A- z)- 1 - (B - zyHdxdy <Ctd. 
k Jc 



□ 



Lemma 3. Let h S S°(e,C) be real-valued and f G A. Then h is a bounded self-adjoint operator 
and 

\\f(h)-f(h)\\=0(e 2 ). (35) 

Proof. The Weyl quantization of a real valued symbol in S°(e,C) is bounded and symmetric and 
thus self-adjoint. Following Exercise 24 of Chapter 2 in [Mar02 , it suffices to show (|3"5"|) for the 
resolvents at z 6 C \ R: 

(h - z)- 1 - V w ((h - z)- 1 ) = (h- z)- 1 ^ -(h-z) Op w ((/i - z)- 1 )) 
= (h-z)- 1 (l-Op w ((h-z)#(h-z)- 1 )) 
= (h - zy'jf Uh - z), (ft - z)' 1 } +Q{e 2 )) . 



□ 
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